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ABSTRACT 

Spectroscopic observations from the Large Binocular Telescope and the Very Large Telescope reveal 
kinematically narrow lines (~ 50 kms -1 ) for a sample of 14 Extreme Emission Line Galaxies (EELGs) 
at redshifts 1.4 < z < 2.3. These measurements imply that the total dynamical masses of these systems 
are low (< 3 x 10 9 Mq). Their large [O Illjsoor equivalent widths (500 — 1100 A) and faint blue 
continuum emission imply young ages of 10 — 100 Myr and stellar masses of 10 8 — 10 9 Mq, confirming 
the presence of a violent starburst. The stellar mass formed in this vigorous starburst phase thus 
represents a large fraction of the total (dynamical) mass, without a significantly massive underlying 
population of older stars. The occurrence of such intense events in shallow potentials strongly suggests 
that supernova-driven winds must be of critical importance in the subsequent evolution of these 
systems. 

Keywords: galaxies: dwarf — galaxies: evolution — galaxies: formation — galaxies: high-redshift 
galaxies: starburst 


1. INTRODUCTION 

The z > 1 universe contains a remarkably large num- 

ber of galaxies with extremely luminous nebular emis- 

sion lines in comparison to their faint blue continua 

(van der Wel et al. 2011). These extreme emission line 

galaxies (EELGs) can have [O III] and/or Ha equiva- 

lent widths (EWs) in excess of 500 A (Atek et al. 2011; 

van der Wel et al. 2011; Brammer et al. 2012b). Such 
observations indicate that young starbursts dominate the 
energy output of these otherwise faint galaxies. While 
similar objects do exist at z < 1 (Cardamone et al. 2009; 
Izotov, Guseva, & Thuan 2011), they have much lower 

comoving number densities and thus their abundance 

must be a strong function of time. 

More massive populations of older stars could easily 
be outshone by the young starbursts: an old stellar pop- 

ulation could have mass-to-light ratios up to ~50 times 

larger than those of the bursts in the near-IR, so the 
main uncertainty in the interpretation of the observa- 
tions hinges on the determination of the total masses of 
these systems. This missing datum prevents us from un- 

derstanding the role of this mode of star formation in the 
broader context of galaxy formation: when these bursts 
occur in truly low-mass galaxies (~ 10 8 Mq), the EELGs 


may represent the main formation mode of present-day 
dwarf galaxies, as argued by van der Wel et al. (2011); 
alternatively, when these bursts occur in more massive 
systems (~ 10 9 Mq) or at higher redshifts, we may be 
witnessing the early formation stage of Milky Way-type 
galaxies. 

Dynamical mass estimates are key in addressing this is- 
sue. For this purpose we now present near-infrared spec- 
troscopy of 14 EELGs at redshifts 1.4 < 2 < 2.3 with 
[O HI] 5007 equivalent widths > 500 A from the Large 
Binocular Telescope and the Very Large Telescope. In 
a companion paper (Maseda et al., in prep.; henceforth 
M13), we describe the observations and data reduction 
in full detail, as well as a more complete interpretation 
of the nature of these objects, while here we focus exclu- 
sively on their mass estimates. Besides the first dynami- 
cal mass measurements of such extreme galaxies, we also 
present improved total stellar mass estimates through 
stringent modeling of the continuum and emission line 
measurements from CANDELS multi-wavelength pho- 
tometry (Grogin et al. 2011; Koekemoer et al. 2011) and 
low-resolution grism spectroscopy from the 3D-HST sur- 
vey (Brammer et al. 2012a). 

We adopt a flat ACDM cosmology with = 0.3 and 
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Ho = 70 km s 1 Mpc 1 throughout. 

2. CANDIDATE SELECTION AND OBSERVATIONS 

A full description of the sample selection, observations, 
and data reduction is given in M13. Here, we focus on 
the 14 objects 15 with restframe equivalent widths > 500 
A in [O HI] 5007 from that emission line-selected sam- 
ple: five are from the photometrically-selected sample 
of van der Wei et al. (2011) in the GOODS-S and UDS 
fields, one is from the grism sample of Straughn et al. 
(2011) in the GOODS-S field, and the remaining ob- 
jects were selected based on their 3D-HST grism spec- 
tra in the COSMOS, GOODS-S, and UDS fields. Al- 
though the targets are very faint in the continuum 
(rriFi 4 o,AB > 24), the emission lines are strong, with 
fluxes > 10 -17 ergs -1 cm -2 A -1 , making emission line 
detections possible with ~1 hour integrations on 8m class 
telescopes. We observe five objects using long-slit ob- 
servations with the X-Shooter wide-band spectrograph 
(Vernet et al. 2011) at the Very Large Telescope from 
August to December 2012 (one slit contained two ob- 
jects), focusing here on the combined Y JHK NIR region 
(1024—2480 nm), although it simultaneously observes in 
the UV-Blue and the Visible regions. Four had 40 minute 
integrations, while one object was observed for a total of 
120 minutes in the near-IR over the course of two nights. 
The remaining ten objects in the sample were observed 
using the LUCI1 multi-object spectrograph (Seifert et al. 
2003) at the Large Binocular Telescope with four sepa- 
rate masks between April 2012 and March 2013 in the 
J-, H-, and/or K - band (depending on the redshift, as 
we targeted [O III] and/or Ha) for a minimum of 45 
minutes per band. Not all targeted objects in the LUCI1 
sample have detected emission lines, either due to their 
intrinsic faintness or contamination from OH skylines: 
in this case, two objects in the total sample had a priori 
equivalent widths greater than 500 A but severe contam- 
ination from OH skylines at the predicted position of the 
lines prevents a line extraction. For all observations, see- 
ing was better than 1" and typically between 0.3" and 
0.8". All exposures were dithered by 3" to decrease de- 
pendence on the pixel-to-pixel detector variations and 
defects. 

Reduction of the X-Shooter data is performed using 
version 2.0.0 of the ESO XSHOOTER pipeline 16 , which 
provides merged, 2D near-IR spectra. Reduction of the 
LUCI1 data is performed using a custom pipeline, with 
the wavelength calibration done using the OH sky lines 
and based on the XIDL routines 17 . For the brightest emis- 
sion lines, we also use XIDL for the final sky subtraction, 
which uses a spline-fitting algorithm to measure and re- 
move the lines. 

Identified emission lines in the ID spectra are fit with 
Gaussian functions, where all lines in a subregion of the 
spectrum (i.e. [O III] AA4959,5007 and H/3) are forced to 

15 There are 15 objects in M13 that fulfill the equivalent width 
criteria. COSMOS- 10320, although fulfilling the criteria, exhibits 
broad and asymmetric [O III] (and also Ha) of 240±13 kms -1 . 
As this object is an obvious outlier (with a potential AGN contri- 
bution), we exclude it from the subsequent analysis and focus on 
the remaining 14 objects. 

16 http : //www. eso . org/ sci/sof tware/pipelines/ xshooter/ 
xsh-pipe-recipes . html 

17 http : //www .ucolick. org/~xavier/IDL/ 
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Figure 1 . Plot of the [O III] 5007 emission line for each object in 
our sample, scaled to the peak flux value. Gray regions show the 
+/- 1-cr flux uncertainties. Typical uncertainties are smaller than 
10“ 4 in redshift and ~8 kms -1 in a. 

have the same width: only the ratio of the two [OIII] 
components is fixed to 2.98 (Storey & Zeippen 2000). 
When both [O III] and Ha are observed for a single ob- 
ject, we take the width of the higher-S/N line complex to 
be the “true” width, which is [O III] for this entire sam- 
ple. The two line widths are always consistent within 
l-cr. 

Extracted emission lines are shown in Figure 1. The 
sample has a median line width of 48 km s -1 with an av- 
erage uncertainty of 8 km s -1 , after correcting for seeing 
and instrumental broadening which is typically < 20% of 
the intrinsic line width. 


3. DYNAMICAL AND STELLAR MASSES 
3.1. Dynamical Mass Measurements 

The velocity dispersions derived above can be used to 
estimate the dynamical masses according to: 

M dyn = C- (1) 

Here, we have adopted the half-light radius r e g as 
the virial radius. We take r e g as the half-light radius 
from van der Wei et al. (2012), who provide size mea- 
surements from the F125W and F160W HST/WFC3 
CANDELS imaging. We choose the filter that does not 
contain the [O III] emission line to ensure that the size 
is measured from the continuum light as much as possi- 
ble. The typical r e g is 1 kpc, which is larger than the 
HWHM of the PSF, so these sources are indeed resolved. 
As noted in Weiner et al. (2006), kinematic estimates 
using line widths yields a variety of results: Rix et al. 
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Table 1 

Summary of Near-IR Observations and Masses 


ID 

RA 

Dec 

Instrument 

Zspec 

5007 

a [OIII] 

( km s -1 ) 

M-dyn 

M* (MAGPHYS) 


(deg) 

(deg) 



(A) 

(Mq) 

(Mq) 

COSMOS-15144 

150.156769 

2.360800 

LUCll 

1.412 

1130T247 

43.3T8.9 

9.11±0.34 

S-lOt^e 

COSMOS-13848 

150.176987 

2.345390 

LUCI1 

1.444 

888T351 

46.7T14.4 

9.22±0.40 

8.58/g; 11 

COSMOS-12807 

150.159546 

2.333301 

LUCll 

1.583 

628T152 

38.2±10.0 

8.88T0.37 

7-95t°; 18 

UDS-7444 

34.473888 

-5.234233 

X- Shooter 

1.621 

713±42 

71.1T5.7 

9.66T0.33 

o yo + 0.07 
° — 0.16 

COSMOS-16207 

150.183090 

2.372948 

LUCll 

1.649 

536±20 

47.7T9.5 

9.40T0.34 

8 43+ 0 - 17 

UDS-3760 

34.428570 

-5.255318 

X- Shooter 

1.664 

731±86 

48.2±5.9 

9.04±0.31 

7 98”*~ 0,11 

' • yo — 0.09 

UDS-3646 

34.426483 

-5.255770 

X-Shooter 

1.687 

701±95 

54.7T6.1 

9.47T0.33 

o n +0.12 
-0.13 

GOODS-S-17892 

53.171936 

-27.759146 

X-Shooter 

1.687 

693±47 

52.3T5.7 

9.05T0.30 

8-95 ±8:1? 

GOODS-S-26816 

53.071293 

-27.705803 

X-Shooter 

1.738 

861±66 

54.4±4.5 a 

8.86T0.31 

8-53t“;°® 

UDS-11484 

34.431400 

-5.212120 

LUCll 

2.185 

723±95 

54.2±9.4 

9.35T0.34 

8.97/p 

COSMOS-11212 

150.124237 

2.313672 

LUCll 

2.199 

598T189 

40.3T8.9 

8.78T0.36 

o 77 +O .23 
' -0.26 

COSMOS-8991 

150.095352 

2.287247 

LUCll 

2.220 

714±85 

30.9T9.0 

8.65T0.40 

9 05”*~°' 21 
U.UO_o 27 

UDS-14655 

34.391373 

-5.195310 

LUCll 

2.297 

503±34 

61.0T10.8 

9.67T0.33 

9 37+ 0 - 11 
-0.31 

UDS-4501 

34.390755 

-5.250803 

LUCll 

2.298 

803T162 

57.8T9.7 

9.07T0.33 

8-32/J].i 9 


Note. — All IDs refer to the CANDELS catalog for that particular field (COSMOS, UDS, or GOODS-S), all equivalent widths are 
quoted in the restframe, and all masses are log quantities. 
a Ho width. 


(1997) calculate C = 2.8 for inclined rotating disks, while 
Barton & van Zee (2001) calculate C = 2.1 for blue com- 
pact dwarfs; Erb et al. (2006) use a simple geometric cor- 
rection to obtain C = 3.4. Here we adopt C = 3, with a 
conservative uncertainty of 33%, as in Rix et al. (1997). 
Note that this value of C would be the same if we as- 
sume that these systems are spherical. We find that the 
14 EELGs have log (Mdy n /M Q ) ranging from 8.7 to 9.7, 
with a median of 9.1 and an average uncertainty of 0.3. 

Sources of systematic uncertainties in these estimates 
affect the results. First, the half-light radius is not neces- 
sarily equal to the virial radius, and some have irregular 
morphologies. Second, these systems likely have a more 
complicated dynamical structure. Third, and most im- 
portantly, the ionized gas in these systems is likely not 
virialized. Dynamical mass estimates on approximately 
kiloparsec scales, which we assume to be dominated by 
the baryonic content of these galaxies, are upper lim- 
its given our assumptions about the dynamics of these 
systems, namely since non-gravitational motion may sig- 
nificantly contribute to the measured velocity dispersion. 

3.2. Stellar Mass Measurements 

With confirmed redshifts, measured EWs of multiple 
lines, and multi-wavelength photometry, we are now in 
a position to estimate the stellar masses and improve 
upon the photometry-only method of van der Wei et al. 
(2011). We take 0.3— 2.2/rm photometry for the six ob- 
jects in the GOODS-S field from Guo et al. (2013) and 
the three objects in the UDS field from Galametz et al. 
(2013): only six of the objects have uncontaminated 
IRAC Ch.1/2 fluxes, so for consistency we do not in- 
clude them for any object, although the mass estimates 
are similar. No such multi-wavelength photometry is as 
of yet available for the six objects in the COSMOS held. 
For these objects we use CANDELS 4-band HST pho- 
tometry (ACS F606W and F814W, WFC3 F125W and 
F160W). 

We follow the procedure of van der Wei et al. (2011) to 
estimate ages and mass-to-light ratios (and hence stellar 


masses) for the EELG sample based on photometrically- 
estimated H/3 EWs and the Starburst99 model (SB99, 
Leitherer er al. 1999). This method assumes a constant 
star formation rate history and equal extinction of the 
emission lines and the continuum. The total EW is 
determined using the excess in F125W (F160W) com- 
pared to the continuum measured from F814W/F775W 
and F160 (F814W/F775W and F125), depending on the 
redshift. The EW is due to the combined effects of [O 
III] AA4959,5007 and H/3; we assume that H/3 contributes 
1/8 of the combined EW. For the output burst age, the 
SB99 model predicts the (restframe) V-band mass-to- 
light ratio, which can be applied to the implied U-band 
luminosity determined from our J- or H-band observa- 
tions, depending on the redshift of the object. Extinc- 
tion is estimated by comparing the UV-continuum slopes 
of the SB99 model and the U/7-band ACS photometry, 
using the Calzetti et al. (2000) extinction law for star- 
burst galaxies. Applying this method with continuous 
star formation, a Chabrier (2003) IMF, a high-mass cut- 
off of 100 Mq, and a metallicity of 0.2 Z & , we derive 
burst masses in the range log (M/Mq) = 7.9 — 9.4 for our 
sample of 14 EELGs. 

With the additional data we have, we can use more so- 
phisticated techniques to estimate stellar masses. Here 
we fit the broadband spectral energy distributions, in- 
cluding line fluxes measured from 3D-HST grism spec- 
troscopy, of our galaxies using a custom version of the 
MAGPHYS code 18 (da Cunha, Chariot, & Elbaz 2008) that 
includes nebular emission computed using the Pacihci et 
al. (2012) model (Chariot & Longhetti 2001, C. Pacihci 
et ah, in prep.). The stellar emission is computed us- 
ing the latest version of the Bruzual & Chariot (2003) 
models using a Chabrier (2003) IMF, and the attenua- 
tion by dust is accounted for using the two-component 
prescription of Chariot & Fall (2000). The nebular emis- 
sion is computed using the CLOUDY photoionization code 
(Ferland 1996), as described in Pacihci et al. (2012). 

18 http : //www . iap.fr/magphys/magphys/MAGPHYS . html 
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Figure 2. Comparison of MAGPHYS- and Starburst99-derived stel- 
lar masses for our sample. Starburst99 utilizes the equivalent 
width of H/3 (determined from photometry alone) to calculate the 
masses, while MAGPHYS utilizes the full photometric SED and the 
emission line fluxes. 

MAGPHYS uses a Bayesian approach to compare the mea- 
sured photometry of observed galaxies with an exten- 
sive library of 100,000 spectral energy distribution mod- 
els spanning a wide range in star formation histories, 
ages, and metallicities. The standard MAGPHYS priors 
(calibrated using more massive galaxies at low redshift) 
are not optimized for this specific population of young 
ages and low metallicities, so we have modified the stan- 
dard priors to include a larger fraction of low metallic- 
ities (between 0.025 and 1 Z©), and younger ages by 
allowing both rising and declining star formation histo- 
ries, all with superimposed random bursts of star forma- 
tion. This method results in stellar masses in the range 
log (M/M©) = 8.0 — 9.4, which are listed in Table 1. 

In Figure 2 we compare the SED-based stellar mass 
estimates with the photometric EW-based stellar mass 
estimates. The photometric method predicts values 1.1 
times (0.05 dex) lower, with an intrinsic scatter (observed 
scatter minus the average uncertainty in the mass) of 0.20 
dex, consistent with no systematic offset. The MAGPHYS 
modeling results reinforce the notion that these galax- 
ies are dominated, in terms of stellar mass, by a very 
young stellar population. While the MAGPHYS modeling 
uses much more information, the crucial elements in both 
mass estimates are the blue continuum and the strong 
emission lines, which strongly constrain any modeling 
approach. Line and continuum luminosities, as well as 
dust content and metallicity, are comparable in the two 
models (see M13). Differences in population synthesis 
models leads to slight intrinsic differences in the M/L for 
a given age and metallicity. In addition, the photometric- 
EW method is very sensitive to the H/3 equivalent width, 
which in this case is simply assumed from the photom- 
etry alone and does not utilize the additional spectro- 
scopic H/3/[0 III] information, also limiting its applica- 
bility to other redshifts. When the ratio is lower than the 



Figure 3. Dynamical masses determined from the velocity width 
of the emission lines versus stellar masses from the MAGPHYS SED fits 
to the full optical/near-IR SEDs for our equivalent width-selected 
sample. The dashed line shows the average value of 27.1% of the to- 
tal dynamical mass made up by stars. The gray point is COSMOS- 
10320 , which is not considered in the analysis. Open diamonds are 
from Erb et al. (2006) for star-forming galaxies at £ ~ 2. Although 
the M dyn values were derived in different manners (see Sec 3.1), the 
relationship between M dyn and M* is similar for the two samples. 

1/8 as used in van der Wei et al. (2011), the resulting 
masses are higher. Finally, the simplistic assumptions of 
constant star formation and fixed metallicity, while ap- 
plicable in general, lead to a very narrow range in ages 
of log(age/yr) = 6.8 — 7.4: the older ages predicted by 
MAGPHYS are coupled with higher masses, where the old- 
est objects have the largest mass differences between the 
two methods. We thus deem the more rigorously-derived 
MAGPHYS masses to be more reliable. 

Figure 3 compares the MAGPHYS stellar mass estimates 
with the dynamical estimates, log (Md yn /M*) = 0.57 
(27% of the total mass is in stars) ± 0.21 (random) ± 
0.30 (sys) for the sample where the 0.30 dex systematic 
uncertainty is from the dynamical mass (see Section 2). 
The 0.21 dex random uncertainty contains the contribu- 
tions from the measurement uncertainties and the limited 
sample size. The three points closest to the Md yn = M* 
line illustrate the challenges to any modeling approach: 
two of them are the only z ~ 2.2 galaxies from the COS- 
MOS sample, where the 4-band CANDELS photometry 
does not sample any continuum redward of [O III], and 
the third is an object with two distinct components in 
the WFC3 imaging, where the assumptions contained in 
the dynamical mass estimate may not accurately reflect 
the true conditions in the system. 

The low dynamical masses confirm the low-mass na- 
ture of these systems directly and exclude the presence 
of large amounts of unseen stars, gas, dust, or dark mat- 
ter that exceed the observed amount of stellar matter 
by more than a factor of five. Our implied maximal gas 
fractions do not exceed those for more massive galax- 
ies at similar redshifts, which range from ~ 30 — 80% 
(Daddi et al. 2010; Tacconi et al. 2013). As seen in Fig- 
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ure 3, our galaxies have similar Md yn /M* ratios to the 
starforming sample of Erb et al. (2006), albeit with EWs 
(and hence specific star formation rates) that are a factor 
of four higher. In M13 we further discuss the implications 
of our constraints under the assumption that the unseen 
matter is in the form of gas or dark matter. 

4. CONCLUDING REMARKS 

In this Letter we show kinematic line widths in the 
range 30 — 70 kms” 1 for a sample of 14 EELGs (with 
EW > 500 A) at redshifts 1.4 < 0 < 2.3. This consti- 
tutes the most comprehensive direct mass measurements 
for such low mass galaxies at these epochs. From this, 
we can now robustly conclude that these galaxies have 
low total masses (~ 10 91 Mq) and that the observed 
young starbursts (~ 10 8 ' 5 M 0 ) constitute a substantial 
part of their total mass budget. These are not high-mass 
(> 10 10 M 0 ) galaxies: they do not contain a signifi- 
cant population of faint, older stars, and they do not 
have extremely high gas fractions. This confirms that 
the abundant population of EELGs aX z > 1 constitute 
a common starburst phase among low-mass galaxies at 
these epochs given the duty cycle arguments presented 
in van der Wei et al. (2011). The wide range in observed 
stellar masses for a given equivalent width could indicate 
that these starbursts do not represent a homogeneous 
population but rather a phase in their evolution, given 
that the highest and lowest mass galaxies in this sample 
are potentially the progenitors of very different galaxies 
today. However, we could also be seeing a homogeneous 
population in which the scatter in mass is driven entirely 
by a scatter in age. A larger data set will help us to dif- 
ferentiate between the two scenarios. 

Given the intensity of the starbursts and the shallow 
potential wells in which they occur, supernova-driven 
winds likely dominate the star formation history and sub- 
sequent evolution of these systems (Larson 1974). The 
starbursts may affect the central dark matter distribu- 
tion (e.g, Navarro, Eke, & Frenk 1996; Read & Gilmore 
2005; Pontzen & Governato 2012; Zolotov et al. 2012) 
and produce cored profiles that are commonly ob- 
served in present-day, low-mass galaxies. For a review 
see de Blok (2010), and Walker & Penarrubia (2011); 
Amorisco & Evans (2012) for recent advances. Our cur- 
rent data set does not allow us to make stronger con- 
clusions about the presence of feedback and winds via 
asymmetric or separate broad/narrow components in in- 
dividual galaxies. However, with future spectroscopic 
studies of these objects, we will be able to search for 
such signals in stacked spectra. 

In the present-day universe, such extreme starbursts 
are very rare (e.g. Cardamone et al. 2009), but at early 
epochs (z > 4 — 6) such events may well be the rule 
rather than the exception. It is becoming increas- 
ingly clear that strong emission lines affect the search 
for and interpretation of high-z galaxies. Strong emis- 
sion lines galaxies at moderate redshifts (z ~ 2) can 
masquerade as drop-out selected z > 10 candidates 
(see discussion in e.g., Coe et al. 2013; Bouwens et al. 
2013; Ellis et al. 2013; Brammer et al. 2013). Further- 
more, for true high-redshift galaxies these strong emis- 
sion lines are likely omnipresent (Smit et al. 2013) and 
affect the broad-band spectral energy distribution, so 
they should therefore be included in the modeling as de- 


scribed here in Section 3.2 (also see Curtis-Lake et al. 
2013; Schaerer, de Barros, & Sklias 2013). However, the 
results presented here are encouraging: we suggest that 
if strong emission lines are evident, then it is likely that 
the total stellar mass does not exceed by much that of the 
young stellar population traced by the blue continuum. 

MVM is a member of the International Max Planck Re- 
search School for Astronomy and Cosmic Physics at the 
University of Heidelberg, IMPRS-HD, Germany. This 
work is based on observations taken by the 3D-HST Trea- 
sury Program and the CANDELS Multi-Cycle Treasury 
Program with the NASA/ESA HST, which is operated 
by the Association of Universities for Research in As- 
tronomy, Inc., under NASA contract NAS5-26555, and 
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